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Abstract 

Fusion cross sections of 28 Si + 28 Si have been measured in a range above the barrier with a very small energy step (AE/ a /, = 0.5 
MeV). Regular oscillations have been observed, best evidenced in the first derivative of the energy-weighted excitation function. 
For the first time, quite different behaviors (the appearance of oscillations and the trend of sub-barrier cross sections) have been 
reproduced within the same theoretical frame, i.e., the coupled-channel model using the shallow M3Y+repulsion potential. The 
calculations suggest that channel couplings play an important role in the appearance of the oscillations, and that the simple relation 
between a peak in the derivative of the energy-weighted cross section and the height of a centrifugal barrier is lost, and so is the 
interpretation of the second derivative of the excitation function as a barrier distribution for this system, at energies above the 
Coulomb barrier. 
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1. Introduction 

Heavy-ion fusion dynamics near and below the Coulomb bar¬ 
rier is a matter of continuing interest, since it allows a deep 
insight into the fundamental problem of quantum tunnelling 
of many-body systems facilitated by channel coupling effects. 
Cross section enhancements, barrier distributions and, more re¬ 
cently, fusion hindrance effects have been observed and are be¬ 
ing investigated. 

Moreover, oscillatory structures were evidenced a long time 
ago in the fusion excitation function of light heavy-ion systems 
like 12 C + 12 C, 12 C + 16 0 and 16 0 + 16 0 ElIlEEl, in the en¬ 
ergy region above the Coulomb barrier. Analogous oscillations 
were found for 20 Ne + 20 Ne 0. In that work, it was suggested 
for the first time that such features are due to successive partial 
waves entering the fusion cross section as their centrifugal bar¬ 
riers are exceeded. In such cases the separation between nearby 
barriers is large with respect to the intrinsic energy width asso¬ 
ciated with the quantal penetration, so that the oscillations may 
become observable. 

This topic was further analyzed in the more recent work of 
Esbensen 0, where earlier data on 28 Si + 28 Si 0 were com¬ 
pared to detailed calculations in the coupled-channels (CC) 
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model, along with the experimental evidences on lighter sys¬ 
tems. Those previous data 0 on 28 Si + 28 Si above the bar¬ 
rier have experimental errors and energy steps just too large to 
allow a clear-cut conclusion about the existence of oscillating 
structures. A more detailed investigation for this system was 
performed by Aguilera at al. in Ref. 0 using the y-ray tech¬ 
nique, where the authors did not observe any oscillation in the 
range E/ a b ^58-100 MeV. 

The existence of such structures, and the detailed effects pro¬ 
duced by channel couplings in the energy dependence of the fu¬ 
sion cross sections, can be best revealed by the first and second 
derivatives of the excitation function multiplied by the energy 
(the so-called “energy-weighted excitation function”) 


£>(£) = 


d(Eo-f) 

dE 


B(E) = 


d 2 (Eo-f) 

dE 2 


( 1 ) 


where E is the center-of-mass energy. As discussed in detail 
in Ref. 0, and following the Hill-Wheeler expression, each 
partial-wave cross section behaves like a Fermi function in the 
no-coupling limit, so that D(E) is the sum of individual cen¬ 
trifugal barriers weighted with the factor (2 L +1), each one 
centered at V&(L). 

The concept of a barrier distribution as defined by B(E) 
has been very helpful in the analysis of several sets of fusion 
data near and below the Coulomb barrier for medium-mass and 
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heavy systems, often giving a ’’fingerprint” of the relevant cou¬ 
pled channels. Its definition was inspired by Wong’s formula 
whose second derivative is a symmetric distribution centered at 
the s-wave Coulomb barrier. 

The barrier distribution B(E) does not carry any information 
about the individual L-dependent barriers. Indeed, when the 
peaks of neighboring L values overlap strongly, one can sum 
up their contributions as in Wong’s formula, and the barrier dis¬ 
tribution is then given by B(E). The overlap between near-by 
peaks, however, diminishes with increasing L, which typically 
results in the breakdown of Wong’s formula above the barrier 
(and of the interpretation of B(E) as a barrier distribution). This 
has been recognized J6l in light, symmetric systems like 16 0 + 
16 0 where only even values of L contribute to fusion and the 
peaks are more distant from each other, but not in the fusion 
of heavy systems, where the condition for separating the indi¬ 
vidual centrifugal barriers (Eq.(ll) of Ref. 0) requires large 
L values where many reaction channels open up and smear out 
the structures. 

The hindrance phenomenon of heavy-ion fusion at deep sub¬ 
barrier energies has been observed for several systems mm 
unna in the last decade. The onset of fusion hindrance has 
often been associated with the energy where the logarithmic 
derivative. 


UE) = 


1 d(Ecrf) 
Eo-f dE 


( 2 ) 


reaches the value (named L cs ) expected for a constant astro- 
physical S -factor Ifl3l . At that energy the S -factor develops a 
maximum as a function of the energy. However in several cases 
the hindrance effect is not strong enough to produce an S -factor 
maximum on. 

The phenomenon is very intriguing and far from being fully 
understood. The structure of the two colliding nuclei El and, 
possibly, couplings to transfer channels mm. affect the en¬ 
ergy threshold below which hindrance shows up. 

The M3Y+repulsion potential of Ref. fl8l produces a rel¬ 
atively shallow potential in the entrance channel and it has 
been capable of explaining the fusion hindrance phenomenon in 
many cases, when applied in CC calculations. There are other 
models on the market that can explain the fusion hindrance phe¬ 
nomenon without employing a shallow potential, for example, 
the model by Ichikawa et al. 03. In that model the hindrance 
is caused by the damping of collective excitations for overlap¬ 
ping nuclei {20). In order to resolve the differences between the 
two models, it is therefore of great interest to test these models 
against new observables, such as the oscillations that appear in 
the measured cross sections at high energy. 

Hindrance effects are more clearly observed in heavier sys¬ 
tems, where, on the other hand, possible consequences of 
stronger channel couplings also deserve attention. An inter¬ 
mediate case like 28 Si + 28 Si calls for interesting investigations. 
High-precision data in a sufficiently wide energy range would 
allow exploiting (and checking) both definitions of Eq.([T]i to 
compare the experimental results with detailed CC calculations. 

The purpose of this work is then twofold: 1) to search for os¬ 
cillations in the fusion excitation function for 28 Si + 28 Si above 


the barrier by careful measurements with very small energy 
steps, and 2) to obtain a consistent interpretation of both the 
sub-barrier fusion excitation function eh and of the oscilla¬ 
tions within the same theoretical CC model. A partial and 
preliminary account of the experimental part of this work was 
given in Ref. [22]]. 


2. Experimental 

28 Si beams with intensities -15-30 pnA were delivered by 
the XTU Tandem accelerator of the Laboratori Nazionali di 
Legnaro of INFN. Targets of 28 Si (with an isotopic enrichment 
of 99.93%) 50 /./g/crrr thick on 15 /rg/cm 2 carbon backings fac¬ 
ing the beam, were used. 

Two separate series of careful measurements of the excita¬ 
tion function (named I and II run) have been performed in the 
energy range ^62-78 MeV (above the barrier), with a step small 
enough (AE/ fl £ = 0.5 MeV) to resolve possible oscillations. For 
each energy, at least 10000 fusion evaporation residues (ER) 
were detected, thus reducing the statistical error to 1% or less. 
The relevance of the accuracy of the 28 Si beam energy in such 
measurement is obvious. A particular care has been devoted 
to this issue, by stepping the field in the 90° analysing mag¬ 
net of the accelerator only downwards, so to minimise possible 
hysteresis effects. The maximum uncertainty in the beam en¬ 
ergy was measured to be +0.13% Il23l (^90 keV at 70 MeV). 
When the energy is monotonically decreased, as in the present 
experiments, the relative beam energy uncertainty is a factor 
3-4 lower. 

The ER were detected at 8i a b- 3° using the same set-up and 
procedures described in Refs. E1E0 , based on an electrostatic 
beam separator. This set-up is very simple to operate, allowing 
fast and reliable measurements of relative and absolute cross 
sections. In the present case, the absolute scale was fixed by 
normalising the relative yields to the cross sections of Ref. Ell 
at corresponding energies. 

Fig. [T shows the excitation function of 28 Si + 28 Si in linear 
and logarithmic energy scales. The excitation function above 
the barrier looks very smooth at first sight, but a closer inspec¬ 
tion reveals small glitches. In the figure, the statistical error bars 
are smaller than the symbol size and the red and blue dots refer 
to the two series of measurements. By extracting the derivative 
of the energy-weighted excitation function |6) D(E) we obtain 
rather regular oscillations as shown in Fig. [2] although the un¬ 
certainties become obviously larger. 

It is worth noting that two well defined peaks at E c m -35 
and 36.5 MeV definitely show up in both series of data (I and II 
run in the figure). This reassures us about the quality of the per¬ 
formed measurements. However, a third peak at E cm ^33.5-34 
MeV is clearly observed only in the II run, due to 2-3 points 
of the I run being significantly lower. In our experience, this 
originates from very small unwanted changes of the beam con¬ 
ditions (focusing, direction) when changing the energy, beyond 
the controls of the diagnostics. The consequences are hardly 
visible in the excitation function (see Fig. [TJa)), but the repre¬ 
sentation in terms of its derivative amplifies the effect. 
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3. Coupled-channels analysis 



Figure 1: (color online) Fusion excitation function of 28 Si + 28 Si in linear (a) 
and logarithmic (b) scales. The results of the CC calculations discussed in the 
text are also shown in (b). 



E (MeV) 
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Figure 2: (color online) The derivative of the energy-weighted fusion excitation 
function D(E) of 28 Si + 28 Si in an energy range above the barrier. The derivative 
is obtained as the incremental ratio between successive points, with an energy 
step AE c/ „ = 0.75 MeV. The lines are the results of the theoretical calculations 
discussed in the text. 


CC calculations have been performed, using the same for¬ 
malism and structure input recently employed for the low- and 
near-barrier data of 28 Si + 28 Si Gil . In the following, as in that 
work, the Chi and ChlO calculations refer to the no-coupling 
limit and to the calculation where 10 coupled channels have 
been considered. The M3Y+repulsion potential m has been 
used, as discussed below. 

The differences from the potential used for the sub-barrier 
data ED are a larger difFuseness of the imaginary potential, a w 
= 0.3 fm instead of a w = 0.2 fm, and a slightly smaller density 
radius of 28 Si, R = 3.125 fm instead of R = 3.135 fm. These 
modifications were made because they improve the fit of the 
ChlO calculation to the data of Ref. ED by reducing the x 2 IN 
from 1.7 to 1.2. The diffuseness of the density that is used to 
construct the M3Y potential, as well as the diffuseness associ¬ 
ated with the repulsion, and the strength of the imaginary po¬ 
tential, have not been varied (a = 0.48 fm a r = 0.398 fm and Wo 
= 5 MeV, respectively). 

The calculated excitation function is compared to the data in 
Fig. |T]b). The corresponding derivative D(E) is shown in Fig. [2] 
(red curve) and is in fairly good agreement with the observed 
oscillations. We point out that the present ChlO calculation 
was calibrated to fit the measured low-energy cross sections of 
Ref. ED- The prediction it makes at higher energies agrees 
well with the new data shown in Fig. [TJb) and the first deriva¬ 
tive of Ecr shown in Fig. [2]is also in fairly good agreement with 
the data. It is remarkable that the whole set of data, including 
the oscillations, are now reproduced within a single theoretical 
model. 

The old data of Ref. Q (not shown here) were originally an¬ 
alyzed by the CC calculations of Ref. j6[ that used a rather large 
radius parameter R = 3.17 fm, and a r = 0.378 fm, with a rela¬ 
tively strong imaginary potential (a w = 0.5 fm. Wo = 10 MeV). 
The blue line in Fig. [2] is that calculation which agrees with 
the present observations rather well. However, that calculation 
does not fit the recent sub-barrier cross sections |2T1 (Fig.jTJb), 
blue curve) at all, even when a weaker imaginary potential is 
used. This leads us to the important conclusion that measuring 
the sub-barrier excitation function is essential for disentangling 
the ambiguities in the choice of the ion-ion potential that arise 
when only considering the data above the barrier. 

The result of the no-coupling calculation Chi using the po¬ 
tential of this work is also reported in Fig. [2] as a black dashed 
line. Weak oscillations can be observed providing a poor fit to 
the experimental data. A strong effect of couplings on the posi¬ 
tion and amplitude of the peaks is obviously seen. It is tempting 
to associate the observed oscillations to the penetration of suc¬ 
cessive centrifugal barriers, as in lighter systems where neigh¬ 
boring barriers are well separated. However, this association 
turns out to be somewhat distorted as discussed below. 

Fig. [3j a) shows again the observed oscillations together with 
the results of the ChlO calculation. The individual contributions 
to D(E) from the angular momenta L = 16-20 are also plotted. 
The comparison shows that the experimental peaks near 33, 35 
and 37 MeV correlate with the peaks produced by the angular 
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Figure 3: (color online) (a) The first derivative of the energy-weighted cross sections. The results of Chi0 calculations are shown. The solid red curve is the full 
result with a maximum angular momentum of L m = 38. The result for L m =14 is also shown and so are the individual contributions for L = 16-20. (b) The results of 
Chi calculations are shown. The experimental peaks at 33, 35, and 37 MeV are associated with L = 14, 16, and 18, while they are primarily associated with L = 16, 
18, and 20 in Chi0 calculations. 


momenta L = 16, 18, and 20. However, the height of these peaks 
constitute less that 50% of the total distribution, partly because 
the peaks overlap and partly because the peaks are fragmented 
as discussed below. It is therefore concluded that it is not pos¬ 
sible to assign a particular angular momentum to each exper¬ 
imental peak. Rather, each peak receives contributions from 
several /.-values. For example, the small peak that appears near 
33 MeV for L = 20 is believed to be a real and not a spurious 
peak because a similar small peak appears for the other values 
of L. However, they are not visible in the linearplot of Fig.[3ja). 

Chi calculations, i. e., without any coupled-channels effects 
(see also Fig. [2]), are reported in Fig. [3j b). The Chi calculation 
does not contain any strong peaks. This is because the peaks 
for L = 14, 16, and 18 in the no-coupling limit are broad and 
overlap, so that their sum is essentially flat. It is also seen that 
the three experimental peaks correlate with the three calculated 
peaks for L = 14, 16, and 18. This correlation, nevertheless, 
does not allow one to assign an angular momentum to the ex¬ 
perimental peaks, because the strong coupled-channels effects 
lower and fragment the effective centrifugal barriers. 

The present data, along with the Chi and ChlO calculations, 
point to the essential contribution of channel couplings for the 
appearance of oscillations. The magnitude of the oscillations is 
quite sensitive to the strength and the diffuseness of the imag¬ 
inary potential. However, the peak positions are insensitive to 
these parameters. The diffuseness a w = 0.3 fm was chosen be¬ 
cause it optimises the fit to the old data of Ref. 11211 . It produces 
oscillations in ChlO calculations that are in reasonable agree¬ 
ment with the new data, which demonstrates a consistency of 
the old and the new data sets. This is valid for 28 Si + 28 Si, 
but it is not a general conclusion. Indeed, in other (even near¬ 
by) systems, coupling effects might as well destroy oscillating 
structures reminiscent of penetration of successive L-barriers. 

We point out that the shallow potential we have used pro¬ 
duces rather thick centrifugal barriers and the associated pene¬ 


tration factor will therefore change quickly as the beam energy 
increases across the barrier height. The rapid rise of the pene¬ 
tration factor results in relatively narrow peaks as illustrated in 

Fig. 0 a). 

4. Fusion barrier distributions 

In general, the barrier distribution one obtains as the second 
derivative of the energy-weighted experimental cross sections 
or of CC calculations (see Eq. (|T])) may contain several peaks, 
due to couplings to other reaction channels during the fusion 
process. The distributions extracted from the Chi and ChlO 
calculations discussed in this paper using the M3Y+rep inter¬ 
action, are compared in Fig.[4ja) to the results of the previous 
measurement SID and Run II of the present measurement. 

The near- and sub-barrier data produce a strongly asymmet¬ 
ric distribution that peaks at 28 MeV. The present data reveal 
a second peak at 32 MeV, in good agreement with the ChlO 
calculation, although the first peak is somewhat over-predicted 
while the second one is well accounted for. 

In this respect it is useful to plot (see Fig. |4jb)) the first 
derivative of the energy-weighted cross section for the angular 
momentum L- 0 which is peaked at the Coulomb barrier in the 
Hill-Wheeler approximation j26j. It resembles the barrier dis¬ 
tribution shown in Fig.jfJ'a), both with respect to the first and the 
second peak obtained in the ChlO calculation. As discussed in 
the Introduction, this is a nice confirmation of the original idea 
of Rowley et al. (25), namely, that the second derivative defined 
in Eq. {T} can be interpreted as a barrier distribution. In partic¬ 
ular, the experimental peak observed at 32 MeV in Fig.[4ja) is 
a real peak of the barrier distribution because it is confirmed by 
both of the ChlO calculations that are shown in Fig. |4ja) and 
Fig-EJb), respectively. 

The experimental barrier distribution shown in Fig.[4|a) de¬ 
velops some structures at energies above the peak at 32 MeV 
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Figure 4: (color online) Barrier distributions defined in Eq. Q and obtained from the present Run II data and from Ref. ed are compared in (a) to distributions 
obtained from the Chi and ChlO calculations discussed in the text. Panel (b) shows the first derivative of the energy weighted cross sections calculated for the 
angular momentum L=0. 


that do not appear convincingly in Fig. |4jb). In that energy 
range, we have evidence that the contributions of successive L- 
values gradually lose overlap. This implies the breakdown of 
Wong’s formula, and, consequently, of the interpretation of the 
second derivative of Eq. 0 as a barrier distribution for a system 
as heavy as 28 Si + 28 Si. 

5. Conclusions 

In conclusion, we have shown the results of detailed mea¬ 
surements of the fusion excitation function of 28 Si + 28 Si 
with very small energy steps, which reveal regular oscilla¬ 
tions, best evidenced in the first derivative of the energy- 
weighted excitation function. We have been able to reproduce 
these high-energy oscillations and the sub-barrier cross sec¬ 
tions, within the same coupled-channels model using the shal¬ 
low M3Y+repulsion potential. 

It appears that the existence of oscillations is tightly bound 
to channel couplings in this relatively heavy system, while in 
lighter cases the oscillations have been suggested to be related 
to the overcoming of successive centrifugal barriers well spaced 
in energy. In 28 Si + 28 Si, the oscillations do appear, but the one- 
to-one relation between each peak and the height of a centrifu¬ 
gal barrier is lost because of strong coupling effects. Checking 
the importance of the oblate deformation of 28 Si in this, calls 
for an analogous experiment on the nearby system 30 Si + 3<) Si 
because 3() Si is essentially a spherical nucleus. 

As the last point of this article we also suggest that, for 28 Si + 
28 Si, the interpretation of the second derivative of the excitation 
function as a barrier distribution breaks down at energies well 
above the Coulomb barrier. 
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